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We investigate the paramagnetic-metal-to-antiferromagnetic-metal and antiferromagnetic-metal-to- 
antiferromagnetic- insulator transitions using a slave-boson mean-field theory. To this effect, we 
discuss the ground state of the half-filled Hubbard model as a function of t' /t and correlation 
strength U, where t and t' are the hopping amplitudes between nearest and next-nearest neighbors, 
respectively. The metal-insulator transition at a critical f/jviIT is of second order for small levels 
of magnetic frustration, t' /t < 0.06, and of first order for large ones, t' /t > 0.06. The insulator 
is always antiferromagnetically ordered, while the metal exhibits a second-order transition from a 
paramagnetic to an antiferromagnetic state up to t' /t — 0.14, as U is increased. We also contrast 
these findings with what we obtain in Hartree-Fock approximation. 
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The correlation-driven metal-insulator transition, or 
Mott transition observed in materials such as V2O3 
[H,D and NiS2-xSc:r H^Uli is a nonperturbative problem 
usually tackled within the Hubbard model of strongly 
correlated electrons. This model describes itinerant elec- 
trons subject to an on-site repulsion U comparable or 
greater than the bare bandwidth 2D. 

In an early work, Brinkman and Rice Q investi- 
gated the Mott transition from the metallic side using 
Gutzwiller's variational scheme In this approxima- 

tion, the metal is described as a strongly renormalized 
Fermi liquid. A low-energy scale ZD [Z is the quasiparti- 
cle residue) collapses linearly in U as the Mott transition, 
occuring at a critical UbKi is approached from the metal- 
lic side. ZD is a measure for the renormalized Fermi 
energy. 

We investigate the implications of antiferromagnetic 
long-range order on either side of the Mott transition. 
To this effect, we introduce magnetic frustration, which 
helps stabilize an antiferromagnetic metallic phase for 
not too large levels of frustration and causes the insu- 
lating side to favor antiferromagnetic long-range order 
P-|l2|. We determine the phase diagram and the or- 
ders of the transitions in mean-field theory using a slave- 
boson technique. Our work complements previous stud- 
ies [p^p^ and adds some analytical insights on the in- 
terplay of electron-electron correlations and magnetism. 

The single-band Hubbard model is given by 
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(1) 



where we take the amplitudes tij to be nonzero only be- 
tween nearest and next-nearest neighbors, in which cases 
they equal t and t' , respectively, and Cio- are cre- 
ation and annihilation operators for an electron of spin 
a at site i, and fiia = cf^Cia- In this work, we consider 
the two-dimensional cubic lattice and restrict ourselves 
to half filling and zero temperature. 



We use a slave-boson method ||T^ in which the approx- 
imation of Gutzwillcr, Brinkman, and Rice is recovered 
on the saddle-point level, but which at the same time 
is open to various generalizations such as the inclusion 
of antiferromagnetic long-range order. The slave-boson 
method yields a ground-state energy in good agreement 
with quantum Monte Carlo simulations including antifer- 
romagnetic order and spiral states , or with exact 
diagonalization studies . It has also been used to de- 
termine a magnetic phase diagram [ p^f and to go beyond 
the Hartree-Fock approximation in problems involving 
complicated spatial structures such as stripes Q . 

In the Kotliar-Ruckenstein approach, two aspects of a 
physical electron are separated: that it is a fermion and 
that it affects the occupancy of some site. The first as- 
pect is taken into account by a fermionic field fi^, while 
the possible occupancies of the sites are described by 
bosonic fields: ei describes empty, pia singly occupied, 
and di doubly occupied sites. The physical electron field 
is represented as Cia- = Ziafi<y with Zi^, = (1 - p^Pjo- - 
dtd^)-'/HetP^^ + Pt.dt){l - e+e, pUp.,^.)-^l\ 
while appropriate constraints eliminate unphysical states 
|l3| . Thus, the problem posed by the Hubbard interac- 
tion is shifted to that of keeping track of the backflow of 
bosonic excitations, z^^Zj^, accompanying the itinerant 
fermions, 

Proceeding along the lines of Ref. we first set 

up the functional-integral representation of the Hubbard 
model in terms of the above-mentioned auxiliary fields, 
integrate out the fermions, and solve the remaining prob- 
lem in the saddle-point approximation. To describe anti- 
ferromagnetism, we divide the lattice into two sublat- 
tices, A and B, and look for solutions satisfying the 
following relations between the sublattice Bose fields: 
CB = BA, Pbg = PA-a, ds = c^A, and m = p\^ - p\^, 
where m is the staggered magnetization. For our re- 
sult, we need the dispersion relations of the renormalized 
quasiparticle bands, 
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= — 4^' COS kx COS ky + rjt^J X'^ + 4(cos k-c + cos ky)^^ 

(2) 



where the lattice spacing has been set equal to one, X 
is some dynamically generated staggered magnetic field, 
and 77 = ±1. The equations for the density per site 
(which at half filling is equal to one) and the staggered 
magnetization, 



(3) 
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N ^ + 4(cos/c^ +cosfcj,)2 

X [/(e,r_m-/i)-/(e,,^[X]-A) 
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can be solved unambiguously for X and the effective 
chemical potential /i, to yield functions /i(m) and X(m). 
In Eqs. (|^) and (|^), N is the total number of lattice sites, 
the sum is over the first Brillouin zone, and /(e) = 6(— e) 
is the Fermi function at zero temperature. From the 
mean-field equation 
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■K{m) + U ^0, 
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where the functions q{m, d?) and K(ra) are given by 
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[1 - + v/(l - 2d2)2 _ ^2] ^ (6) 
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we obtain the average portion of doubly occupied sites as 
a function of the staggered magnetization, <Plrn). This 
function along with Eqs. (H) and (|^) allow to write the 
ground-state energy per site as a function of the staggered 
magnetization as 



e(m) — q{m)K{m) + U(f'{m). 



(8) 



This result has an intuitive interpretation: K{m) is 
the kinetic energy of noninteracting lattice fermions with 
nearest and next-nearest neighbor hopping, subject to 
an internal staggered magnetic field tX(m). The renor- 
malization factor q{m) accounts for the reduction of the 
hopping amplitudes due to the local correlations and is 
characteristic of the Gutzwiller approximation. In our 
scheme, q arises from the expectation value (z^Zja-) and 
thus represents the average effect of the backflowing slave 
bosons. The second term in Eq. (^ is the contribution 
of the Hubbard interaction to the energy. 

For J7 = 0, Eqs. (|) and (|) yield d'^ = (1 - m'^)/A and 
q = 1, so e(TO) = K(m). 



For = m = 0, Eqs. (|)-(0) imply: X = jl = 0; K = 
2 de Z?o(e)e, where Do{e) is the density of states for 
noninteracting electrons; = 1(1 — jj-) with Uc — Si\K\; 
and q = I ~ (y-)"^ . We thus recover a result of Ref. [ p^ . 

In the strong-correlation limit, U t, we find up to 
leading order in t/U: d"^ = At^/U'^, q = I - At^/U'^, 
K = ~8t'^/U, and TO = 1 - 8t'^/U^. Hence, the ground- 
state energy is e = —At^/U. This result is qualita- 
tively correct: In the strong-correlation limit, the half- 
filled Hubbard model can be mapped onto the Heisen- 
berg antiferromagnet with an exchange coupling con- 
stant J = At'^/U. If we treat the electron spins as Ising 
spins and use that the coordination number of our two- 
dimensional lattice is four, we also obtain e = — J. Ther- 
fore, the mean-field theory gives rise to the energy scale 
J and correctly accounts for a ground-state energy of the 
order of — J per site. Finally, we note that the leading 
1/U corrections do not depend on t'. More generally, 
the entire insulating phase is unaffected by next-nearest- 
neighbor hops, t' enters only via Eqs. (^), (^ and (|^). 
For the insulator, however, only the band — but the 

entire one — contributes to the k sums in these equa- 
tions, and t' drops out. 

The ground state for given model parameters U and 
a = t' /t corresponds to the minimum of the energy 
function (H). a is a measure for the degree of mag- 
netic frustration and is varied from zero to one. We 
may restrict ourselves to positive a, since a sign change 
of a is tantamount to a particle-hole transformation, 
Qirj Gxp{iQRi}a^_^ (a and a"*" describe the quasiparti- 

cles with dispersion (||), Q = (tt, tt), and Ri is the lattice 
vector to site i). The evolution of e(TO) as a function of 
U and a reveals how the transitions between the various 
phases take place. 



For a 7^ 0, The inverse susceptibility, x 



IS 



known explicitly and changes its sign from posi- 
tive to negative as U is increased to above Uq = 
8\Km{Y\K{0)\/{\K{0)\-t^) - 1). Here, I/7 = 
{l/N)J2k^/\coskx -|-cosfcj,| and J^k denotes the sum 
over those regions of the first Brillouin zone which are 
restricted by the condition ef:_{X = 0) < //(to = 0) < 
(^^^{X — 0). Uq and 7 depend only on a. 

Whether the system is metallic or insulating depends 
on the value of the ground-state magnetization: If it ex- 
ceeds a certain value, TOmit, a gap opens up in the single- 
particle spectrum and the system goes insulating. This 
can be seen from Eq. (^ if we use that to increases mono- 
tonically as a function of X . Consequently, the insulator 
is always antiferromagnetically ordered. We infer from 
Eqs. (||)-(^ that tomit does not depend on U, but on 
the level of magnetic frustration: Due to perfect nesting, 
limc^o ™mit = 0. As a is turned on, tomit increases 
monotonically as a function of a. 

In our numerical investigation, we used a tetrahedron 
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method [T^] and, for a < 0.2 and m < mMiT, up to 10^ x 
10'* points to do the k sums in Eqs. (|), (|) and (0). In the 
discussion of our numerical results, we must distinguish 
between three regimes of magnetic frustration and may 
restrict the discussion of e(m) to positive magnetizations. 
Fig. |l] illustrates how e{m) evolves as a function of U and 
a. The resulting phase diagram is displayed in Fig. |^. 

For small levels of magnetic frustration, < a < 0.06 
(first column of Fig. |l|) , we first find a second-order tran- 
sition from the paramagnetic to the antiferromagnetic 
metal at a critical value J7mag = Uq. Upon further in- 
creasing [/, the resulting minimum is continuously shifted 
towards higher magnetizations until it crosses mMiT at 
a second critical value, J7mit- Consequently, the metal- 
insulator transition is also of second order. Since tomit 
vanishes as a — > 0, both transitions coincide in this limit, 
C^mag = C^MiT — Uq- Wc also kuow that Uq = if a = 0. 

For intermediate levels of magnetic frustration, 0.06 < 
a < 0.14 (middle column of Fig. the magnetic transi- 
tion is still of second order and occurs at C/q, t^mag = Uq. 
On the other hand, ttimit is sufficiently large for the 
metal-insulator transition to take place differently: Be- 
fore the ground-state magnetization of the antiferromag- 
netic metal reaches mMiT upon increasing [/, a second 
minimum at a magnetization above ttt-mit has emerged 
and become the absolute minimum of e(m). Conse- 
quently, the metal-insulator transition is now of first or- 
der. The transition lines Umit{cx) from the small-a and 
intermediate-a regimes meet at a = 0.06 (filled circle 
in Fig. H). At this point, the two degenerate minima of 
the first-order transition merge at mMiT- Finally, the 
antiferromagnetic metallic phase disappears gradually as 
a 0.14 (filled square in Fig. ||). 

For large degrees of magnetic frustration, a > 0.14 
(right column of Fig. |^), the antiferromagnetic and metal- 
insulator transitions coincide, Umit = C^mag, because the 
second-order transition at Uq is now preempted by the 
first-order one: By the time the minimum at m = 
bifurcates, the one above tomit has already evolved into 
the absolute one, and remains to be so, as U is further 
increased. 

We have also considered the Hubbard model in the 
Hartree-Fock approximation. Earlier work on the t — t' 
Hubbard model in Hartree-Fock approximation was car- 
ried out in Rcfs. 0, |ll|, and 1^. As for the two- 
dimensional half-filled case, these works did not conclu- 
sively answer whether an antiferromagnetic metal is sta- 
ble in a certain parameter regime. 

Within our formalism, the Hartree-Fock approxima- 
tion turns out to be tantamount to taking q = 1 and 
d?(m) — (1 — m^)/4 in Eq. (||), while the function 
K{m) is again determined by Eqs. (||), (||), and (|^). 
As a consequence, = 2(7i — f//4), which changes 
its sign at Uo{a) = At^{a). This is the small-a limit 
of the corresponding mean-field expression. Expanding 



e^^(m) one step further yields the same — |mplog|m| 
term as in mean- field theory, but with q = 1. This sug- 
gests that the magnetic transition is of second order and 
takes place at UQ{a). In fact, our numerical investiga- 
tion of the evolution of e^^(m) as U is increased re- 
veals the following behavior for all finite values of a: Af- 
ter a second-order transition from the paramagnetic to 
the antiferromagnetic metal at f/mag = Uq{q), the min- 
imum of e^^(TO) is rapidly but continuously displaced 
towards higher magnetizations, until it exceeds ttt-mit- 
Thus, the metal-insulator-transition is also of second or- 
der. Fig. H displays the phase diagram in Hartree-Fock 
approximation. The error bars of the transition lines 
^^mag(a) and UMTT{oi) become of the same order than 
UMwiot) — t^mag(a) for a < 0.001, and our numerics is 
trustworthy down to about a ~ 0.002. 

In summary, we have investigated the effect of mag- 
netic frustration on the metal-insulator transition in 
the two-dimensional half-filled Hubbard model within a 
slave-boson approach, and we have compared our results 
to the Hartree-Fock approximation. Within the slave- 
boson mean-field theory, our main results are: First, 
magnetic frustration helps stabilize an antiferromagnetic 
metal for t' /t < 0.14. Second, for t' /t > 0.06, the metal- 
insulator transition is of first order. Finally, all other 
transitions between the various phases are of second or- 
der. By contrast, in Hartree-Fock approximation, the 
magnetic and metal-insulator transitions are always sep- 
arate and of second order. Both in Hartree-Fock approx- 
imation and in the slave-boson mean-field theory, the in- 
sulator is always antiferromagnetically ordered. 

We gratefully acknowledge discussions with R. Fresard 
and R. Chitra. This work was supported by the NSF 
DMR 95-29138. E.L. is partly funded by the Deutsche 
Forschungsgemeinschaft . 
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FIG. 1. The function e(m) as a function of U and the mag- 
netic frustration. The columns from left to right correspond 
to the small-, intermediate-, and large-a regimes, respectively. 
The respective values are a = 0.02, 0.1, and 0.5, and corre- 
spond to JTiMiT ~ 0.09, 0.26, and 0.61, respectively. Each col- 
umn displays how e(m) changes qualitatively upon increasing 
U, from the paramagnetic metal (top row) to the antiferro- 
magnetic insulator (bottom row). The middle row shows ex- 
amples in the antiferromagnetic metal (first two plots), and 
one after the metal-insulator transition has taken place but 
before the local minimum at m = turns over into a local 
maximum (plot on the right). The plot in the center has 
its minimum at a nonzero magnetization, which is not dis- 
cernible. 




a 



FIG. 2. The various phases and transitions as a function 
of a = t' /t and U. The small-17 phase is the paramag- 
netic metal, while the large-[/ phase is the antiferromag- 
netic insulator. In between, an antiferromagnetically ordered 
metallic phase is sandwiched that disappears at a tricritical 
point (filled square). The dotted and filled lines indicate sec- 
ond- and first-order transitions, respectively. The filled circle 
marks the point where the metal-insulator transition changes 
its order. 
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FIG. 3. The phase diagram in Hartree-Fock approxima- 
tion: The paramagnetic metal (small-f/ regime) and antifer- 
romagnetic insulator (large-f/ regime) are separated by a tiny 
range of U values corresponding to an antiferromagnetically 
ordered metal. Both transitions are of second order. 
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